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ABSTRACT 

We present an analysis of the mechanics of thin streams, which are formed following the tidal 
disruption of cold, low-mass clusters in the potential of a massive host galaxy. The analysis 
makes extensive use of action-angle variables, in which the physics of stream formation and 
evolution is expressed in a particularly simple form. We demonstrate the formation of streams 
by considering examples in both spherical and flattened potentials, and we find that the action- 
space structures formed in each take on a consistent and characteristic shape. We demonstrate 
that tidal streams formed in realistic galaxy potentials are poorly represented by single orbits, 
contrary to what is often assumed. We further demonstrate that attempting to constrain the 
parameters of the Galactic potential by fitting orbits to such streams can lead to significant 
systematic error. However, we show that it is possible to predict accurately the track of streams 
from simple models of the action-space distribution of the disrupted cluster. 

Key words: methods: analytical - methods: numerical - Galaxy: kinematics and dynamics - 
Galaxy: halo - Galaxy: structure - galaxies: interactions 



1 INTRODUCTION 

Recent years have seen a tremendous advance in the quality and 
quantity of observational data for substructure in the halo of our 
Galaxy. Of particular note is the outstanding success of the Sloan 
Digital Sky Survey (SDSS, York et al. 2000), which has uncov- 
ered large numbers of streams in the halo of the Milky Way 
(Odenkirchen et al. 2003; Majewski et al. 2003; Yanny et al. 2003; 
Belokurov et al. 2006, 2007; Grillmair 2006; Grillmair & Dionatos 
2006; Grillmair & Johnson 2006; Grillmair 2009; Newberg et al. 
2009). These streams appear to be the relics of tidally stripped glob- 
ular clusters and dwarf galaxies that have fallen into the Galaxy, 
although it is often the case that no obvious progenitor object can 
be associated with a particular stream. 

It has been empirically noted that such tidal streams appear to 
delineate the orbit of the progenitor object from which they formed 
(McGlynn 1990; Johnston et al. 1996). A consequence of this as- 
sumption is that measurements of separate segments of a single 
stream can be equated to sampling different phases of a single or- 
bit. If either complete phase-space information (Willett et al. 2009), 
or a certain subset of it (Jin & Lynden-Bell 2007; Binney 2008; 
Eyre & Binney 2009a), is known for a sequence of phases along 
a single orbit, it is possible to diagnose the Galactic potential with 
exquisite precision (Binney 2008). Thus the assumption that such 
streams delineate orbits promises tremendous diagnostic power and 
had been frequently invoked (Binney 2008; Eyre & Binney 2009a; 
Odenkirchen et al. 2009; Willett et al. 2009; Koposov et al. 2010; 
Eyre 2010a). 
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Such techniques are hampered by the considerable difficulty in 
obtaining observational data of adequate precision for stars in dis- 
tant streams. More worrisome, however, is the mounting evidence 
that streams do not precisely delineate orbits (Dehnen et al. 2004; 
Choi et al. 2007; Montuori et al. 2007; Eyre & Binney 2009b)— in 
these circumstances, which we will demonstrate to be generic — 
any conclusions drawn from the fitting of orbits to streams could 
be systematically in error. 

It is therefore of vital importance for the continued application 
of such techniques that we obtain an understanding of the relation- 
ship between tidal streams and the orbits of the stars that comprise 
them. In particular, we must determine under what circumstances 
tidal streams delineate orbits, by what measure they are in error 
when they do not, and what can be done to correct this error. There 
has not to date been an exposition of stream formation that funda- 
mentally addresses these issues. 

Most studies have either focused on N-body simulations (Choi 
et al. 2007; Montuori et al. 2007), the confusion of which makes the 
predominant physics hard to isolate, or have attempted to describe 
the problem in terms of conventional phase-space coordinates and 
classical integrals (Dehnen et al. 2004; Choi et al. 2007), which 
makes the problem intractably hard. The work of Choi et al. (2007) 
made some progress towards understanding the dynamical struc- 
ture of clusters at the point of disruption, and they provide a quali- 
tative picture of the evolution of tidal tails, understood in terms of 
classical integrals. However, they are unable to make predictions 
for stream tracks on the basis of this picture alone, and they are 
ultimately forced to rely on N-body simulation for quantitative re- 
sults. In this paper, we approach the problem using action-angle 
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variables, in which the physics of stream formation finds a natural 
and simple expression (Helmi & White 1999; Tremaine 1999). 

We confine our investigation to the formation of long, cold 
streams, such as may form from tidally disrupted globular clusters. 
We do so for two reasons. First, a low mass for the progenitor clus- 
ter simplifies the description of its orbital mechanics, because of 
the lack of dynamical friction and other feedback effects in its in- 
teraction with the host galaxy. Second, thin, long streams provide 
the strongest constraints upon the Galactic potential, because any 
orbit delineated by them can be observationally identified with less 
ambiguity (Binney 2008; Eyre & Binney 2009b). Hence, it is long, 
cold streams that are of primary interest for use as probes of the 
potential. 

We study the mechanics of stream formation immediately fol- 
lowing the tidal disruption of a progenitor cluster. In much of 
the work that follows, the assumption is made that stream stars 
feel only the potential of the Galaxy, i.e. the stream stars do not 
self-gravitate. This assumption is generally a fair one: the stars in 
streams are generally spaced too widely for their self-gravity to be 
of consequence (Dehnen et al. 2004). Indeed, we will demonstrate 
below that self-gravity becomes negligible shortly after stars are 
stripped from our model clusters. 

The remainder of the paper is arranged as follows: Section 2 
discusses the action-angle variables in which we perform our anal- 
ysis. Section 3 discusses the basic mechanics of stream formation 
and propagation. Section 4 discusses the detail of stream forma- 
tion in spherical systems, with the aid of some examples. Section 5 
relates the action-angle structure of streams to the corresponding 
real-space manifestation. Section 6 examines the action-space dis- 
tribution of disrupted clusters using N-body simulation. Section 7 
describes the consequences of optimising potential parameters un- 
der the faulty assumption that streams follow orbits. Section 8 dis- 
cusses stream formation in flattened systems, using an axisymmet- 
ric Stackel potential as an example. Finally, Section 9 presents our 
concluding remarks. 



2 ACTION-ANGLE VARIABLES 

We use action-angle variables to analyse stream formation and 
propagation. The usefulness and theoretical basis of action-angle 
variables is extensively discussed in §3.5 of Binney & Tremaine 
(2008; hereafter BT08). 

Action-angle variables are a set of canonical coordinates, like 
Cartesian phase-space coordinates, that can be used to describe sys- 
tems in Hamiltonian mechanics. Actions J are useful because they 
constants of motion. Moreover, the angle variables that are con- 
jugate to them evolve linearly in time 



0{t) = 9(0) + fit, 



where we have introduced the frequencies 
dH 



n i (j) = e i = —. 



(i) 



(2) 



Hence, the motion of a system described by action-angle variables 
is very easy to predict. 

§3.5 of Binney & Tremaine 2008 discusses at length the cal- 
culation of action-angle coordinates in various systems. Here, we 
note that standard methods for calculating action-angle coordinates 
require the Hamilton- Jacobi equation to separate (§3.5.1 of BT08). 
This condition is met by all spherically symmetric potentials, but 
excludes any non-spherical potential that is not of Stackel form (see 



Section 8.1 below). Given separability, the action corresponding to 
a coordinate q is given by 



l - j p g dq, 



(3) 



where the integral is over the closed path that encloses a single 
oscillation of the coordinate q along an orbit. 

We note that just as spherical systems are naturally described 
by spherical polar coordinates (r, 0), the natural actions for such 
systems are ( J r , L), where J r is the radial action, and L — J$ + 
\L Z \ is the angular momentum. Throughout this paper, and without 
loss of generality, when spherical systems are under study we set 
L — L z and confine the motion to the (x, y) plane. 

In the sections that follow, we utilize the equations from 
§3.5.2 of BT08 to transform between the action-angle coordinates 
( J r , L, 9 r , 6^) and conventional phase-space coordinates, when in- 
vestigating spherical potentials. The action-angle coordinates used 
when investigating non-spherical potentials are discussed in Sec- 
tion 8. 1 below. 



3 THE FORMATION OF STREAMS IN ACTION-ANGLE 
SPACE 

Consider a low-mass cluster that has recently been tidally dis- 
rupted, with just enough time having passed so that the stripped 
stars no longer feel the effects of one another's gravity. The cluster 
is on a regular orbit, identified by its actions J () , in a fixed back- 
ground potential, which has a Hamiltonian H(3). Suppose further 
that in the locality of Jo, the Hamiltonian is well described by the 
Taylor expansion, 



H(3) = H + fi ■ <5J + i,5J T D S3, 
where S3 = 3 — Jo, D is the Hessian of H 
D * H 



and ft is the frequency vector of the cluster's orbit 



dH 



The frequency vector ft of a nearby orbit J is then 

n(j) = n + D-<5J. 



(4) 



(5) 



(6) 



(7) 



If the disrupted cluster has some spread in actions AJ and angles 
AOo, corresponding to the spread in position and velocity of its 
constituent stars, then the spread in angle-space after some time t 
is given by (Helmi & White 1999; Tremaine 1999; BT08 §8.3.1), 



A0(t) = tACl + A9 ~ tASl, 



(8) 



where the near equality is valid when t Aft 3> A9o, and where we 
have introduced the spread in frequencies, Afi, which are related 
to the spread in actions via 



Aft = D ■ AJ. 



(9) 



In the absence of self-gravity, the action-space distribution AJ is 
frozen for all time. D and ft are functions of J only, and so are sim- 
ilarly frozen. The secular evolution of a disrupted cluster is there- 
fore to spread out in angle-space, with its eventual shape deter- 
mined by Aft = D • AJ, and its size growing linearly with t. 

For a given AJ, what does Aft look like? Equation (8) and 
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equation (9) show that D acts as a linear map between a star's posi- 
tion in action-space and its position in angle-space. D is a symmet- 
ric matrix, so it is characterised by its mutually orthogonal eigen- 
vectors e n (n — 1,3) and its real eigenvalues A n . 

Consider a cluster whose stars are distributed isotropically 
within a unit sphere in action-space centred on J . The correspond- 
ing angle-space structure, resulting from the mapping of this sphere 
by D, will be an ellipsoid with semi-axes of length t\ n and direc- 
tion e„. 

If the An are finite and approximately equal, such an isotropic 
cluster will spread out in angle-space with no preferred direction. 
Eventually, the cluster will uniformly populate the whole of angle- 
space. In real-space the cluster will uniformly fill the entire volume 
occupied by the orbit Jo: such a cluster would not form a stream. 

If one of the A n is much smaller than the others, then D will 
map the cluster into a highly flattened ellipsoid in angle-space. In 
real-space, the cluster will occupy a 2-dimensional subspace of the 
orbital volume of Jo. The precise configuration of this subspace is 
likely to be complex, but it would not form a stream. 

If two of the A„ are small, then D maps the cluster into a line 
in angle-space. In this case, the resulting real-space structure will 
be a filament. The density of this filament will fall linearly with t. 
In a real galaxy, such a structure may therefore persist with a sig- 
nificant overdensity for some time. It is this case that describes the 
formation of tidal streams (§8.3.1 of BT08) investigated in detail in 
this paper. Finally, if all the \„ are zero, there will be no spread, 
and even an unbound cluster will remain extant indefinitely. 

We note that there is no a priori reason for any of the \ n to be 
small. The existence of D imposes no conditions on H in general, 
save that it must be twice differentiable near to Jo. We can write a 
Hamiltonian for which, for a particular J at least, the A„ and the 
e„ take any specified values. It must therefore be a peculiar prop- 
erty of realistic galactic potentials that causes disrupted clusters to 
form streams. 



3.1 Validity 

Let us consider briefly the validity of the preceding analysis. The 
key assumption is that the Taylor expansion (4) holds, which it does 
if 



d 3 H 



-SJk- 



(10) 



3 dJk 3 OJidJjdJk 

If H is dominated by some low power of J, this condition becomes 

SJi < J % . (11) 

In general then, we expect our analysis to be valid if the spread in 
action of the stars in the cluster is small compared to the actions 
themselves, which is likely to be true for cold clusters on high- 
energy orbits around massive hosts. The condition (10) is met in 
detail for all the examples considered in this paper. 



3.2 The problem 

We have seen that the condition for a stream to form is that one of 
the eigenvalues A n of the linear map, D, must be much larger than 
the other two. Herein, we will number the A„ and their correspond- 
ing e„ such that 



Ai > A2 ^ A3. 



(12) 



Under what conditions does the direction of this stream point along 
the progenitor's orbit? 



Taken over a short period of time — say, the time taken for an 
individual star to travel along the length of a stream — every star 
in the stream will follow an almost identical trajectory. If phase is 
disregarded and the trajectories superimposed, they would overlay 
almost perfectly. Streams form because slight differences between 
the trajectories of individual stars grow secularly over the lifetime 
of a stream, which may be many orbital periods of the individual 
stars. There is no fundamental dynamical reason why the track of 
the resulting stream should delineate the progenitor's orbit, nor the 
orbit of any one of the stars, nor any orbit in the governing potential. 

In action-angle coordinates, the trajectory in angle-space of 
an orbit J is given by the frequency vector, fl(J). A necessary 
condition for a stream to delineate precisely the progenitor orbit is 
that the long axis of the angle-space distribution A0 is aligned with 
the progenitor frequency flo. However, we shall see in Section 5.3 
that because real-space position x(J, 0) is a function of J as well 
as angle 0, this condition alone is not sufficient to ensure real-space 
alignment between streams and orbits. 

We note that the angle-space distribution, A0, depends on 
both the potential, D, and the action-space distribution, via AJ, 
and in practical cases it is necessary to consider both when deter- 
mining its gross alignment. 

The problem of stream formation and evolution is patently a 
complicated one. To further our analysis, we will examine sepa- 
rately the role of the linear map D, the real-space map x(J,#), 
and the action-space distribution A J. Having done so, we will pro- 
ceed to use this understanding to examine the formation of streams 
in practical cases. 



4 THE LINEAR MAP D IN SPHERICAL CASES 

In order to analyze the behaviour of the linear map, let us restrict 
ourselves to a case in which our cluster is both isotropic and small 
in A J. The angle-space distribution A0 is now an ellipsoid with 
semi-axes of length A„, and with the semi-major axis of the distri- 
bution aligned with the principal direction of the linear map, ei. In 
this scenario, the stream will be delineated by a single orbit if the 
condition 



ei = fl , 



(13) 



is met, since AJ is small and so x(J, 0) ~ x(Jo, 0). 

We will proceed by examining the forms of fl and D, and 
hence the angle-space geometry of streams, for some example 
spherical potentials. We remind ourselves that such systems can 
always be completely described by two actions. 

4.1 Kepler potential 

There are remarkably few potentials of interest for which the 
Hamiltonian can be written as a function of J in closed form. One 
such potential is the Kepler potential 



$ w = — — > 

for which the Hamiltonian is, 

(GM) 2 



(14) 



(15) 



^ (J) ^ 2(J„ + L)' 

where J r is the radial action, and L is the angular momentum. The 
frequencies can be worked out by direct differentiation, giving 



fl 



(GM) 2 
2(J r + L) V 1 



(16) 
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Table 1. Actions and apses for selected orbits in the spherical potentials 
used in this paper. We generally choose = L so all orbits remain in the 
(x, y) plane. 
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Figure 1. The solid line shows the orbit Kl (Table 1) in a Kepler potential, 
on which a cluster of 50 test particles has been evolved. The particles were 
released at apocentre. The red dots show the positions of the test particles 
near apocentre, after 24 complete orbits, at t = 4.02 Gyr. The blue dots 
show the same test particles near pericentre, approximately half an orbit 
later. In both cases, the dots delineate the progenitor orbit precisely. 



The linear map D is calculated via a further round of differentia- 
tion, and has eigenvalue/vector pairs 



[Ai.Si] = 
[A 2 ,e 2 ] = 



Ai, 












Aa, 













-2(GM) 

(Jr + L) 
-1 



0. 



(17) 
(18) 



The linear map D has only one finite eigenvalue, and so the Kepler 
potential will always form streams from disrupted clusters. Further, 
we observe that since ei oc 17 in all circumstances, streams will 
always perfectly delineate progenitor orbits in Kepler potentials, 
regardless of the action-space distribution of the cluster from which 
they form. In this potential streams exhibit secular spread strictly 
along the orbit, and do not grow wider with time. 

4.1.1 Numerical tests 

We have confirmed the above prediction numerically. We use a Ke- 
pler potential with mass M = 1.75 x 10 M@, which reproduces 
a fiducial circular velocity v c = 240 km s _1 at the approximate so- 
lar radius of i?o = 8 kpc. A cluster of 50 test particles was created, 
sampled from a Gaussian distribution in action-angle space, defined 
by a j — lkpckms -1 and ag — 5 x 10 radians. This cluster 
was placed on the orbit Kl, given in Table 1, in the aforementioned 
Kepler potential. The orbit is highly eccentric, with a pericentre 
radius of about 1.5 kpc and an apocentre of about 13 kpc. 

The cluster was released at apocentre, and evolved for 
4.02 Gyr, equal to 24 complete orbits, by integrating the equations 



Figure 2. Angle-space configuration for the particles shown at apocentre 
in Fig. 1. The solid line shows the frequency vector f2o, with which the 
stream particles are perfectly aligned. Notably, there is no secular spread in 
the direction perpendicular to the stream motion, as predicted in Section 4. 1 . 



of motion for each particle. Fig. 1 shows the real-space configura- 
tion of the particles at the end of this time, and Fig. 2 shows the 
configuration of the same particles in angle-space. 

The cluster has elongated to form a stream that Fig. 1 shows 
to cover approximately half the orbit when the centroid is at peri- 
centre. The stream does not spread in width in either real-space or 
angle-space. In both figures, the stream delineates the cluster's orbit 
perfectly. Fig. 1 shows that this is true irrespective of the real-space 
location of the centroid, so the prediction of the preceding section 
is validated. 

4.2 Spherical harmonic oscillator 

The spherical harmonic oscillator potential 

*(r) = \n 2 r\ (19) 

where Q. is a constant, applies for motion within a sphere of uniform 
density (§3. la of BT08), and is therefore of relevance to galaxy 
cores in the absence of a black hole. The Hamiltonian is 



H(J r ,L) = Q{L + 2J r 
and the frequencies are 



(20) 



(21) 



while D is a null matrix, implying that the eigenvalues are identi- 
cally zero. 

From this latter fact, we conclude that clusters in harmonic 
potentials will always remain in the same configuration in angle- 
space, and will not spread out. Hence, streams cannot form in har- 
monic potentials. We further note that, in any case, it would be 
difficult to tidally strip a cluster in a harmonic potential, since the 
tidal force divide across a cluster 
<9 2 $ 



dF tld 



Or- 



-dr 



(22) 



is independent of galactocentric radius r. Thus, a cluster that is 
bound at apocentre in such a potential will remain bound elsewhere 
along its orbit. 

4.3 Isochrone potential 

The isochrone potential (§2. 2. 2d of BT08) is a simple potential 
which has several useful properties. It behaves as a harmonic oscil- 
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Figure 3. Details of the stream geometry in the isochrone potential of Sec- 
tion 4.3. Upper panel: the misalignment angle in degrees, between the 
principal direction of D and CIq, shown as a contour plot against the ac- 
tions of the progenitor orbit. In all cases, 1? is between 1.2 dcg and 2.2 dcg 
in angle-space. Lower panel: The ratio of the eigenvalues A1/A2. The ra- 
tio is > 10 everywhere and rises sharply with increasing L. The actions 
shown in both plots cover a range of interesting orbits, which are described 
in Table 2. 



lator in the limit of small radius, and as a Kepler potential at large 
radius, thus providing a reasonable model for a spherical galaxy 
across all radii. The form of the potential is 

GM 



$(r) = - 

V ; b + Vb 2 + r 2 

where b is a constant scale. The Hamiltonian is 

(GM) 2 



H(J) 



2[J r + ±(L + VL 2 + 4GMb)] 2 



(23) 



(24) 



To proceed, we require the frequencies, which are obtained by par 
tial differentiation of equation (24) 

(GM) 2 



{Jr + ±(L + VL 2 +4GMb)}^ 



1 + 



(25) 



(26) 



VL 2 +AGMb / 

To construct D, we also require the derivatives of the frequencies 
with respect to the actions, which are obtained directly from equa- 
tions (25) and (26). 

The algebraic forms of both fl and D are untidy, and little 
progress can be made by simply inspecting them. Instead, we pro- 
ceed by working numerically with a specific example. We use the 
isochrone potential with the parameters M = 2.852 x 10 11 M©, 
6 = 3.64 kpc, which is chosen to have a rotation curve maximized 
with v c — 240 km s _1 at the fiducial solar radius of Ro = 8 kpc. 

What then is the geometry of streams formed in this potential? 



Table 2. The coordinate extrema of selected orbits from Fig. 3, illustrating 
the variety of orbits covered by that figure. The actions are expressed in 
kpc km s~ 1 , while the apses are in kpc. 
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Fig. 3 shows the misalignment angle between fio and ei, given by 

■d — arccos (tto ■ ei^j , (27) 

as well as the ratio of the eigenvalues A1/A2, both as functions of 
J. The range of J shown in Fig. 3 covers a variety of interesting 
orbits, described in Table 2. 

The upper panel shows that the principal direction of D is mis- 
aligned with the progenitor orbit by 1-2 deg for all values of J. The 
misalignment is at a minimum for both low energy and high energy 
circular orbits, and at a maximum for eccentric orbits with a guid- 
ing centre close to r = b. The lower panel shows that the ratio of 
the eigenvalues A1/A2 varies from 15 to 35 across the range, with 
the ratio maximized for high energy circular orbits, and minimized 
for high energy plunging orbits. 

On this basis, we expect an isotropic cluster of test particles 
in this potential to form a stream in angle-space that is misaligned 
with fio by 1-2 deg and is 15-35 times longer than it is wide. 
Notably, such a stream is precluded from being delineated by any 
single orbit, close to Jo, in the correct isochrone potential. 

4.3.1 Numerical tests 

We created a cluster of 150 test particles, randomly sampled from 
a Gaussian distribution in action-angle space, defined by aj = 
0.2 kpc km s -1 and a s — 10~ 3 radians. This cluster was placed 
on the orbit II from Table 1, which has an apocentre radius of 
~ 13 kpc and a pericentre radius of ~ 5 kpc. 

The cluster was released at apocentre, and evolved for 94 com- 
plete azimuthal circulations, equal to a period of t = 22.75 Gyr. 
Fig. 4 shows the angle-space configuration of the particles after 
this time. The arrowed, black line shows the orbit of the underlying 
cluster, S7o- The dashed line shows a straight line fit to the distribu- 
tion of particles, which is clearly misaligned with the black line. We 
further note that, unlike the stream in the Kepler potential shown in 
Fig. 2, this stream has clearly increased in width since inception. 

We can predict the shape of this distribution precisely. Plotted 
as a red ellipse in Fig. 4 is the angle-space image under D of a cir- 
cle in action-space. After some time, the angle-space distribution of 
an isotropic cluster of test particles should take the form of a scaled 
version of this image. We see that the image and the particle distri- 
bution are indeed comparable, and that the dashed line is perfectly 
aligned with the principal axis of the image. 

How does this misalignment manifest itself in real-space? 
Fig. 5 shows the real-space configuration of the cluster at the end 
of the simulation. The orbit of the cluster is drawn with a solid 
black line, while the dashed line is a quadratic curve that has been 
least-squares fitted to the particle data. 

The progenitor orbit is clearly a poor representation of the 
stream. Although the orbit passes through the centroid of the 
stream, as expected, the curvature of the orbit is too low to match 
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Figure 4. The angle-space distribution of a cluster of test particles, evolved 
on orbit II in the isochrone potential of Section 4.3.1 for 94 complete az- 
imuthal circulations. The particles are shown at apocentre. The frequency 
vector f2o of the progenitor orbit is shown with an arrowed black line. The 
stream is slightly misaligned with f2o; the dashed line is a straight line fit 
to the positions of the test particles, and clearly demonstrates this misalign- 
ment. Also plotted with a red solid line is the image in angle-space of a cir- 
cle in action-space, mapped by D. The shape and orientation of the image 
reflects the A n and e n of D for this orbit. The ellipse is clearly misaligned 
with the underlying orbit, but is perfectly aligned with the stream particles. 
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Figure 5. Real-space configuration for the stream of test particles shown in 
Fig. 4, plotted in polar coordinates. The arrowed black line shows the tra- 
jectory of the progenitor orbit. The stream formed by the dots falls away in 
radius faster than does the orbit, in both forwards and backwards directions: 
the stream clearly does not follow the progenitor orbit. The dashed line is 
a quadratic curve least-squares fitted to the stream, which shows that the 
stream has a substantially greater curvature than does the underlying orbit. 
The solid blue line is the track predicted by mapping the dashed line from 
Fig. 4 into real-space: it agrees perfectly with the stream. 



the stream adequately. Thus, the small misalignment in angle-space 
is manifest as a significant change in stream curvature at apocentre. 

Since we know the orientation of the long axis of the angle- 
space distribution for this stream, we can predict the stream track 
in real-space. The blue line in Fig. 5 shows the mapping into real- 
space of the dashed line from Fig. 4. The mapping is done by solv- 
ing numerically for the real-space roots of the equations that re- 
late action-angle variables to real-space coordinates: for spherical 
potentials, the appropriate equations are given in §3.5.2 of BT08. 
Fig. 5 shows that, unlike the progenitor orbit, this mapped line de- 
lineates the stream perfectly. 

Although we have only considered the unrealistic case of a 
stream formed from an isotropic cluster, we have nonetheless pro- 



duced a stream which is misaligned from its progenitor orbit in both 
angle-space and real-space. 

We conclude that, in the isochrone potential, the track of the 
stream makes a poor proxy for the orbit of its stars, and that in 
general, streams in this potential cannot be relied upon to delineate 
orbits. 



5 MAPPING STREAMS FROM ACTION-ANGLE SPACE 
TO REAL SPACE 

The example of Fig. 5 shows that an angle-space misalignment can 
alter the curvature of the real-space stream. In this section, we ex- 
amine why this is so, and show the range of real-space effects that 
are a consequence of angle-space misalignment. We will also ex- 
amine the effects of a finite action-space distribution on the real- 
space track of a stream. 



5.1 Non-isotropic clusters 

Up until this point, we have considered our streams to form from a 
cluster of particles that is isotropic in J, resulting in a stream that is 
perfectly aligned in angle-space with the principal direction of D. 

It is not obvious that this is a fair assumption. The structure 
in angle-space, given by equations (8) and (9), is linearly depen- 
dent upon the action-space distribution that generates it. By prop- 
erly choosing that distribution, we can create streams of arbitrary 
shape in angle-space. Nature does not create clusters with arbi- 
trary action-space distributions, so arbitrary-shaped streams do not 
emerge. But what kind of action-space distribution should be con- 
sidered reasonable? 

In Section 6 below, we will investigate the action-space dis- 
tribution of real clusters, by means of N-body simulation. Here, 
we only require a qualitative understanding, in order to know what 
kind of gross action-angle structures we should leam to map. Since 
the action-space distribution fundamentally arises from the random 
motion of stars within the cluster, it is unlikely to be dominated by 
complex structure. We will therefore assume that the gross struc- 
ture is ellipsoidal. But what should be the axis ratio of this ellipse? 

Consider a cluster with isotropic velocity dispersion, a, that is 
on an orbit with apocentre r a and pericentre r p , where it is tidally 
disrupted. In our spherical system, the radial action is given by the 
closed integral 



Jr ~ 27V 



p r dr, 



(28) 



where the integration path is one complete radial oscillation. Now 
consider a particle whose radial momentum p r differs from that of 
the cluster average by Sp r ~ a. We can take a finite difference over 
equation (28) and thus obtain an expression for the the difference 
in radial action between the particle and the cluster 



SJ r ~ — 5p r Ar ~ — aAr, 

TT TT 



(29) 



where Ar 



is the amplitude of the radial oscillation. 



Now consider another particle, whose azimuthal velocity differs 
from that of the cluster by Svt ~ o. The difference in angular 
momentum between this particle and the cluster is 



SL ~ r p 6vt 



(30) 



where we have performed our calculation at pericentre, because 
that is where the cluster is stripped. For the purposes of this section, 
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Figure 6. The upper panel shows a selection of ellipses in action-space, 
each of axis ratio 10, but oriented in different directions. The lower panel 
shows the image in angle-space that results from mapping each of the 
action-space ellipses with the linear map D, calculated for the orbit II in the 
isochrone potential of Section 4.3.1. In the lower panel, the arrowed black 
line is the frequency vector f2o; in the upper panel, the arrowed black line is 
the inverse map of the frequency vector, D _1 f2o. We see that regardless of 
the shape in action-space, the mapped images are all elongated and roughly 
aligned with the orbit, although the alignment is generally not perfect. In 
this example, the misalignment of each of the red and green images is about 
10 dee. 



we are interested in the relative size of the spread in radial action 
AJ r and the spread in angular momentum AL for a disrupting 
cluster. We see that their quotient 



AJ r 
AL 



Ar 

7rr D 



(3D 



Although this ratio will take on every value between (0, oo) as we 
move from a circular orbit to a plunging one, for the orbits likely to 
be occupied by stream-forming clusters, it will typically be within 
an order-of-magnitude of unity. 

The upper panel of Fig. 6 shows a set of ellipses, with axis ra- 
tio 10, placed at various orientations in action-space: an axis ratio 
of 10 would not be considered untypical by the arguments of the 
previous paragraph. The lower panel of Fig. 6 shows the images 
under D of these ellipses in angle-space when evaluated on orbit 
II in the isochrone potential of the previous section. All the images 
are both elongated and roughly oriented towards the principal di- 
rection. We conclude that the images of most action-space ellipses 
under this map — and thus, most streams formed in this potential — 
would be highly elongated and oriented to within a few degrees 
of the principal direction, which is itself oriented to within a few 
degrees of the frequency vector flrj- 

Since the ratio of the eigenvalues for this orbit is ~ 17, it is 
not possible to produce an image in angle-space that is not elon- 



gated towards the principal direction, by mapping an action-space 
ellipse of axis-ratio 10. We note from Fig. 3 that, for this poten- 
tial, the ratio of eigenvalues does not vary much, and nor does the 
principal direction stray from flo by more than a few degrees. We 
therefore conclude that reasonable action-space distributions will 
always result in the formation of streams in this potential, and that 
such streams will always be oriented in angle-space to within a few 
degrees of f2rj- 



5.2 The mapping of action-angle space to real space 

The upper-left panel of Fig. 7 shows three trajectories in angle 
space. The black line is the trajectory of II in the isochrone po- 
tential given in Section 4.3.1. The red line has a frequency ratio 
Q. r /Q.^, that is 10 per cent lower than the black line, and is there- 
fore rotated from it by ~ 2.9 deg. Conversely, the blue line has a 
frequency ratio that is 10 per cent higher than the black line, and is 
therefore rotated from it by ~ 2.5 deg. 

The red and blue lines were chosen to represent likely streams 
in angle-space that could form in the isochrone potential, given the 
results of the previous section. We note that the red line has retarded 
radial phase (relative to the black line) on the leading tail, and ad- 
vanced radial phase on the trailing tail. Conversely, the blue line 
has advanced radial phase in the leading part, and retarded radial 
phase in the trailing part. 

How do these lines map into real-space? The upper-right panel 
of Fig. 7 shows the real-space curve obtained from the lines in the 
upper-left panel of that figure, having chosen the point of inter- 
section such that the lines are phase-matched near apocentre. The 
curves in the upper-right panel have been drawn by assuming that 
all points along each line have the same J. Thus, this plot represents 
the real-space curves of streams oriented in angle-space accord- 
ing to the upper-left panel, but formed from clusters of vanishingly 
small A J. 

In the upper-right panel, we see that the red line has system- 
atically lower curvature than the black line. Conversely, the blue 
line has systematically greater curvature than the black line. This is 
because the red curve has retarded radial phase on the leading tail, 
and advanced radial phase on the trailing tail, and thus is flattened 
with respect to the orbit. Similarly, the blue line has advanced radial 
phase on the leading tail, and retarded radial phase on the trailing 
tail, and thus appears curved with respect to the orbit. 

The lower-left panel of Fig. 7 shows the same lines, but now 
phase-matched at pericentre. Similarly to the upper-right panel, the 
red line again appears less curved with respect to the orbit, and 
the blue line appears more curved with respect to the orbit. The 
lower-right panel of that figure also shows the same lines, phase- 
matched at a point well away from apsis. In this case, a misalign- 
ment between the stream and f2o in angle-space manifests itself as 
a real-space directional misalignment, rather than a curvature er- 
ror as at apsis. This occurs because, unlike at apsis, the mapping 
between angle-space and real-space plane polar coordinates is rel- 
atively undistorted near this point. 

In summary, realistic angle-space distributions in an isochrone 
potential, which may be misaligned with the progenitor frequency 
by a few degrees, will produce real-space streams that exhibit dif- 
fering curvature from the progenitor orbit when observed at apsis, 
and differing directional alignment from the orbit when observed 
far away from apsis. 
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Figure 7. The upper-left panel shows three trajectories in angle space. The black line is the trajectory of orbit II in the isochrone potential of Section 4.3.1. 
The red line has a frequency ratio Q r /Q,f, that is 10 per cent lower than II. The blue line has a frequency ratio that is 10 per cent higher than II. The lines 
are phase-matched at their point of intersection. We note that the red line has retarded radial phase (relative to the black line) on the leading tail, and advanced 
radial phase on the trailing tail. Conversely, the blue line has advanced radial phase in the leading part, and retarded radial phase in the trailing part. The other 
three panels show the real-space images of these three lines, when phase-matched: (upper-right) near apocentre, (lower-left) near pericentre, and (lower-right) 
at a point far from apsis. In each of these three panels, the centroid of the potential is marked with a a cross. 



5.3 Are trajectories insensitive to small changes in J? 

To this point, all real-space tracks have been derived from streams 
in angle-space under the assumption that at all points along that 
stream the action is that of the progenitor, Jo. This assumption is 
only strictly valid in the case of a vanishingly small action-space 
distribution, and for asymptotically large time since disruption of 
the cluster. If a mapping into real-space from a line in angle-space 
is made under this assumption, then a stream generated by a suffi- 
ciently broad action-space distribution will not be accurately rep- 
resented, even though the representation in angle-space may be ex- 
act. This is because the small changes in action that give rise to the 
small changes in frequency also cause small changes in real-space 
trajectory as well. 

When computing a stream track in real-space, it is possible 
to correct for this effect. By inverting equation (9) and eliminating 
Att using equation (8), we find that for a star separated from a 
fiducial point on the stream by angle 80, the difference in action 
between the star and the fiducial point, 5J is given by 



S3 = — D _1 <50, 



(32) 



where td is the time since the star and the fiducial point were co- 
incident. We may therefore guess the correction S 3 for a star's true 
action (Jo + <5J) from its position in the stream, provided we know 

We typically take the fiducial point to be the centroid of the 
stream, following which we may assume td to be the time since 



the first pericentre passage of the cluster on its present orbit. Al- 
though this assumption neglects the possibility that the star could 
have been torn away during a subsequent pericentre passage, we 
note that during tidal disruption, it is the fastest moving stars which 
become unbound. The cluster core that remains after a pericentre 
passage therefore has lower velocity dispersion. Stars subsequently 
torn from that cluster will therefore have a smaller distribution 
in action-space. Consequently, the stars with the largest 86 from 
the centroid — i.e. those for which the 83 correction will be most 
important — must have been torn away at the earliest time, and so 
the assumption that td equals the time since the first pericentre pas- 
sage remains good. 

Just how important is this effect? For small changes in J r , the 
trajectory changes we discuss are expressed as changes in the radial 
amplitude Ar, while the guiding-centre radius r g , which is purely 
a function of L, is held constant. We can estimate the magnitude of 
the effect as follows. 

Consider a cluster on an orbit that is close to circular, whose 
radial action is given by equation (28). Orbital energy E is con- 
served, so close to apsis r = ro, the radial momentum p r is given 
by 



E=~p 2 r (r)+$ eS (r) 

- ^Pr{r) + $ e ff (To) + (r - r ) —j— 



(33) 
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where we have defined the effective potential 

$ eS {r) = $(r) + L 2 /2r 2 . 
Since at apsis p r — 0, 
E = $ e ff(r ). 

Hence from equation (33) we have 



(34) 
(35) 



Pr(r) = J-2(r-r ] 



dr 



= V2(r-r )F cB (r ), (36) 



where we have defined the effective force, 
F eS {r') = -d$> eS /dr\ r ,. 

If (r a , r p ) are apocentre and pericentre respectively, then 



Prif) 



x/2|Feff(r a )|(r a -r-) if r~r a 



v /2|F e ff('r p )|(r - r p ) ifr~r p . 
Hence, we define a global approximation to p r 



Pr(r) 



V2\F oS \(r-r p )(r*-r) 
v/Ar 



(37) 



(38) 



(39) 



where Ar — r a — r p and F e ff a constant set equal to the value of 
F ei r (r) at either apsis, since we assume it takes approximately the 
same value at both. We note that 



/ V( r - r p)( r a - r) dr = -Ar 2 



(40) 



so when combined with equation (39), equation (28) for J r yields 

(41) 



Ar 3 



Jr = lV2\Fci 
O 

We can deduce the value of F e g as follows. We note that 



$ H 



(42) 



and that its derivative is 

d<f>eff _ d<E> L 2 
dr 



dr dr r 
If the rotation curve is relatively flat, then F e ti(r) evaluated at r a 

r g + Ar/2 is 

__vl L 2 
r 



(43) 



Feff : 



d$eff 



dr 



= -A 



Ar 

2 3JJ 
2 + or 4 



(r g + |Ar)3 
_ L^_ f 1 _3Ar 
rl\ 2 r g 
_5w|Ar 
2L 2 ' 



(44) 



and similarly when evaluated at r p , but with opposite sign. Equa- 
tion (41) then becomes 



V5v 2 Ar 2 
Jr ~ SL ■ 

Differentiating the above expression, we find 
d.J r V5v 2 Ar 



(45) 



dAr 



4L 



V5JrVc 

2L 



(46) 



Hence, for a small change 8J r we can estimate the corresponding 
change in the radial amplitude, 5Ar, which is likely to be a good 



estimate for the positional error we would make in assuming that a 
star with action J r + SJ r actually had action J r . 

We can confirm the predictions of the above equations numer- 
ically. The left panel of Fig. 8 shows the real-space trajectory of the 
circular orbit 12, with radius r = 8 kpc, in the isochrone potential 
of Section 4.3.1. Also plotted is the trajectory of an orbit that has 
identical L to 12, but J r = 0.21 kms -1 kpc. 

Clearly, equation (46) ceases to have meaning when faced 
with orbits very close to circular, so we rely for our estimate on 
the integral form equation (45) instead. Equation (45) predicts 
Ar = 0.16 kpc for this perturbation from circular, which appears 
from the left panel of Fig. 8 to be a reasonable estimate. 

From equation (46) we see that the magnitude of the effect 
diminishes with increasing J r as SAr ~ 1/Ar ~ The 
centre panel of Fig. 8 shows this to be the case. The panel shows 
two trajectories in the isochrone potential: one for the orbit 13, and 
one for the same orbit with J r incremented by 0.1 per cent. Equa- 
tion (46) predicts a change <5Ar ~ 3 pc. Close inspection of the 
trajectories confirms an actual SAr ~ 3pc, so the prediction is 
correct, but as is clear from Fig. 8, corrections of such magnitude 
are negligible. 

What follows is an estimate for the positional error made by 
incorrectly guessing L. Consider again a cluster on an orbit close 
to circular. The angular momentum of the cluster is related to the 
guiding-centre radius r g and the circular velocity v c by 

L = i> c r g . (47) 

Consider now a star whose angular momentum is suddenly reduced 
by SL. This star is now at apocentre, since its guiding-centre radius 
has been reduced by 

SL 



Sr g = 



(48) 



where we have assumed that the rotation curve is flat. The pericen- 
tre radius will have been reduced by of order twice the change in 
guiding-centre radius, hence we may write 

28L 



SAr = 



(49) 



Thus, we can predict the change in radial amplitude SAr for a small 
change in angular momentum SL, which is likely to be a good esti- 
mate for the positional discrepancy we would encounter in assum- 
ing that a star with angular momentum L + SL actually had angular 
momentum L. 

Again, we can check the predictions of this expression numer- 
ically. The right panel of Fig. 8 shows the trajectories of two orbits 
in our isochrone potential. The red line is orbit 12, while the black 
line is the same orbit, but with the angular momentum increased by 
1 per cent. Equation (49) predicts Ar ~ 0.16 kpc for this change. 
Fig. 8 shows that this estimate is close to exact. 

The only example we have considered thus far for which this 
effect could be of consequence is the isochrone-potential stream 
shown in Fig. 5. For this cluster, equation (46) predicts that a posi- 
tional error of ~ 2 pc would be accrued by assuming all stars have 
the same radial action. Similarly, equation (49) predicts that a po- 
sitional error of ~ 1.5 pc will be accrued by assuming that all stars 
have the same angular momentum. These errors are insignificant, 
so no corrections are required in this case. 

However, the scale of the AJ distribution in the example of 
Fig. 5 was deliberately chosen to be very small, with a velocity 
dispersion a ~ 4 x 10~ kms -1 . The equations (46) and (49) 
show the magnitude of the trajectory anomaly rises linearly with 
the scale of A J. For a large cluster with a ~ 20kms _1 on the 
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Figure 8. The effect of small changes in actions on real-space orbital trajectories. The left panel shows that a small change in J r (6J r = 0.21 km s — 1 kpc) 
to a circular orbit (12 in Table 1 shown in red) produces a change of <5Ar ~ 0.2 kpc in radial amplitude. The centre panel shows the effect of the same 
perturbation 5J r on an eccentric orbit (13 in Table 1) with the same angular momentum as the circular orbit, but with J r = 207 kms -1 kpc. The effect 
on radial amplitude is so small as to be invisible in this plot. The right panel shows the effect of an increment in angular momentum L of 1 per cent on the 
trajectory of the orbit 12. In this case, the radial amplitude changes by SAr ~ 0.15 kpc. 



orbit 12, this would imply trajectory errors of order ~ 1 kpc, which 
are not negligible. Therefore, consideration of the effects of a finite 
action-space scale on the real-space track of streams is necessary 
in practical cases. 

In the case of a large cluster, even a stream that is perfectly 
aligned with its frequency vector in angle-space will not be delin- 
eated by the progenitor orbit in real-space. The correction described 
by equation (32), to account for variation in action down the stream 
will be required to correctly predict the stream-track from its angle- 
space distribution. In general we would not know the time since 
first pericentre td accurately, although given Jq and a knowledge 
of the extent of the stream, we could make a reasonable guess as to 
its value. However, even a poor, but finite, guess for the value of td 
would likely produce a more accurate real-space stream track than 
would assuming J = Jo everywhere along the stream. 



6 THE ACTION-SPACE DISTRIBUTION OF DISRUPTED 
CLUSTERS 

To this point, we have relied upon the qualitative estimate from Sec- 
tion 5. 1 for what the action-space distribution of a disrupted cluster 
might be. In this section, we investigate the action-space distribu- 
tion of various cluster models using N-body simulation. We further 
utilize our N-body models to confirm the misalignment between 
streams and orbits, and to demonstrate that we can accurately pre- 
dict the real-space track of the stream. 

The action-angle coordinates, of the host-galaxy potential, 
have limited usefulness when applied to the particles in a bound 
cluster, because the actions are not constant with time. Nonetheless, 
they remain a valid set of canonical coordinates and can be legiti- 
mately used to describe the phase-space distribution of the cluster. 
Moreover, we shall see that a disrupted cluster gives rise to a char- 
acteristic distribution in action-space, from which the track of the 
stream can be predicted. 

Fig. 9 shows a segment of each of the orbits 14 and 15 in the 
isochrone potential of Section 4.3.1. These orbits were chosen to 
be fairly representative of those occupied by tidal streams in our 
Galaxy: they have apocentre radius ~ 20 kpc and are moderately 
eccentric to allow for efficient tidal stripping. For our investigation, 
we wish to launch model clusters on each of these orbits, where 
the otherwise-stable cluster has been chosen such that its outermost 
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Figure 9. Plan views of the orbits used in this section. The top panel shows 
14, with (a), (b) and (c) marking the positions of the cluster corresponding to 
the upper-left, upper-middle and upper-right panels of Fig. 10 respectively. 
The bottom panel shows 15, with (a) marking the position of the cluster cor- 
responding to the bottom-right panel of Fig. 12. In both panels, the potential 
in use was the isochrone potential described in Section 4.3.1. 

stars will be torn away by tidal forces when close to pericentre. The 
process by which we choose our model clusters is detailed in the 
next section. 

6.1 Cluster models 

We work with King models (King 1966, BT08) for our clusters, 
since these simple models are both easy to generate and are fairly 
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Figure 10. The action-space distribution of particles for the N-body cluster model C 1 , at different times along the orbit 14. From left-to-right and top-to-bottom, 
these times are: shortly after release; first pericentre passage; first apocentre passage; second apocentre passage; 7th apocentre passage; 14th apocentre passage. 
The solid black line is the inverse map of the frequency vector, D _1 f2o - The dashed line in the bottom-right panel represents a least-squares linear fit to the 
particle distribution. 



representative of some observed globular cluster profiles (Fig. 6. 18, 
Binney & Merrifield 1998). The profile of a King model can be 
defined by W = Qq/<j 2 , the ratio of the central potential to the 
squared- velocity parameter a 2 . For a given W, all resulting models 
are similar in terms of p(f/fo)/ po, where po is the central density 
and fo is the core radius. An exact model is specified by choosing 
po and fo, in addition to W, either directly or through a relation 
with another parameter. 

For a given orbit, we specify our models as follows. Following 
the argument of Dehnen et al. (2004), we note that a cluster of mass 
M c orbiting at a galactocentric radius r from the centre of a host 
galaxy with circular velocity v c , will be tidally pruned to the cluster 
radius fade, where 



We freely choose a profile parameter W and a cluster mass M c , and 
we also specify a galactocentric stripping radius r s > r p , where r p 
is the pericentre radius of the orbit concerned. We then set ft, the 
cluster truncation radius, equal to f t jj c from equation (50), where 
r — > r s and v c — > v c (r s ). The resulting cluster will remain intact 
while r ^> r s , but will have its outermost stars tidally stripped 
when r ~ r p . 

6.2 The disruption of a cluster 

The low-mass cluster model C 1 (Table 3) was specified for the orbit 
14 (Table 1) according to the schema in Section 6.1. We chose a 
low value W = 2 for the profile parameter of for our basic cluster 



model, in order to ensure the presence of many particles near the 
cluster truncation radius ft during successive stripping events. 

A 10 4 particle realization of the cluster model CI was made 
by random sampling of the King model distribution function (equa- 
tion 4. 1 10 of BT08). This cluster was placed at a point shortly after 
apocentre on the orbit 14 in the isochrone potential of Section 4.3.1. 
The cluster was evolved forward in time in the aforementioned po- 
tential by the FVFPS tree code of Londrillo et al. (2003), using a 
time step of dt — td yn /100 and a softening length e as specified 
in Table 3. The simulated time period was 4.81 Gyr, or almost 14 
complete radial orbits. 

Fig. 10 shows the evolution of the action-space distribution 
of the cluster model CI as a function of time. In all the panels of 
that figure, an arrowed black line shows the mapping of the fre- 
quency vector from angle space into action space, D _1 f2o. This 
vector indicates the direction that maps onto Sla in angle-space, so 
any action-space distribution that is aligned with this vector will be 
aligned with tto in angle-space. 

In all cases involving the isochrone potential, D _1 f2o is ori- 
ented exactly along the J r axis. We can understand this from equa- 
tion (26), which shows that in the isochrone potential, the frequency 
direction f2o is a function of L only, and is independent of J r . Thus, 
a line of constant f2o must map into action-space as a line of con- 
stant L. This is a peculiar feature of the isochrone potential, and is 
not true for a general potential. 

The upper-left panel of Fig. 10 shows the configuration of the 
cluster immediately after release, at position (a) in the upper panel 
of Fig. 9. The distribution is ellipsoidal without additional sub- 
structure, which we expect since the distribution in action results 
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Table 3. Details of the cluster models used in this section. The defining parameters for each model are (W, M c ,r s ), while (r'y , cr, t^ yn , e) are derived 
parameters. 
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entirely from the approximately spheroidal density profile, and the 
approximately isotropic velocity dispersion, of the cluster. We un- 
derstand the ellipticity of the action-space distribution from Sec- 
tion 5.1. Indeed, the prediction of equation (31) of A J r / AL ~ 0.3 
for this orbit can be seen to be approximately correct. 

In the upper-middle panel of Fig. 10, the cluster has now 
moved from its point of release to its first pericentre passage, 
marked as position (b) in Fig. 9. We see that the ellipse has flat- 
tened somewhat, and has rotated anticlockwise. In the upper-right 
panel, the cluster has now progressed to the subsequent apocentre 
passage, marked as position (c) in Fig. 9. Here, the cluster is again 
flattened, but now it has rotated clockwise. 

We can qualitatively understand this behaviour as follows. 
Consider a particle in a cluster near apsis. What changes to the ac- 
tions will be made by perturbations to the velocity of this particle? 
A perturbation Sv to the transverse velocity will cause a change to 
the angular momentum 

8L = rSv. (51) 

By means of the mechanism described by equation (49), this SL 
will cause a change in the guiding-centre radius r g , which will 
cause a corresponding change in the radial action, according to 
equation (46). Conversely, a perturbation to the radial velocity will 
cause negligible change to the radial action, since 

8E ~ p r 5p r = rSv ~ 0, (52) 

and J r (E,L) remains unchanged. Hence, the distribution in both 
J r and L is governed primarily by the transverse velocity when the 
cluster is at apsis, and the high degree of correlation observed in 
Fig. 10 reflects this. 

We can confirm this analysis by examining Fig. 11. The upper 
panel shows (v r ,v t ) for the cluster near its first pericentre pas- 
sage, while the lower panel shows the same for the cluster near 
its subsequent apocentre passage. The particles coloured red are 
those with L < 3110 kpckm s _1 and those coloured blue have 
L > 3140 kpckms . The boundary between the colours is very 
sharp in the v t direction, as one would expect, since the particles 
have been coloured on the basis of L. However, we note that par- 
ticles with the full range of v r contribute to both the blue and red 
regions equally, despite Fig. 10 showing that J r is very different for 
these two regions. Hence J r must be independent of v r near apsis. 

We complete our explanation of the orientation of the action- 
space distribution by taking note of the sign of the change in J r near 
apsis. Increasing L will always increase r g , and when the cluster is 
at pericentre, this pushes r g further away, and so increases the ra- 
dial action. Conversely, when the cluster is at apocentre, increasing 
r g brings it closer to the cluster, and so decreases the radial ac- 
tion. Thus, we see that near pericentre, particles with high L will 
have high J r and the distribution will be rotated to have a positive 
gradient A J r / AL. Conversely, at apocentre, particles with high L 
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Figure 11. Radial velocity versus tangential velocity for the CI cluster 
model, on the 14 orbit, at (top panel) the first pericentre passage, and (bot- 
tom panel) the first apocentre passage. 

will have low J r , and so the distribution will be rotated to have a 
negative gradient A J r / AL. 

We can predict the value of the gradient of this distribution, by 
combining equation (46) and equation (49). We find 

where the sign of the radical depends on the apsis under consider- 
ation, as detailed above. We can see from the upper-right panel of 
Fig. 10 that this equation predicts approximately the correct gradi- 
ent of dJ r /dL ~ —0.5, when evaluated for the orbit 14. We further 
note that this equation implies that the gradient will be steeper for 
an orbit of greater eccentricity. 

Let us again examine Fig. 10, and note that the black particles 
are those still bound to the cluster, while the red particles are those 
that are unbound, where we have defined 'bound' to mean all parti- 
cles that are within r t ide of the cluster barycentre. In the upper-right 
plot, which shows the cluster configuration at the first apocentre 
passage, the few unbound particles form an approximately horizon- 
tal distribution. These are particles that have been stripped from the 
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cluster near to the preceding pericentre passage, which is depicted 
in the upper-middle panel of Fig. 10. 

We note two things. First, the bulk of the stripped particles 
have large A J from the cluster centroid; this is simply a conse- 
quence of the high-speed stars being most likely to be stripped. 
Second, we note that although the red particles in the upper-right 
panel span approximately the same range in L as the black parti- 
cles in the upper-middle panel, they span a range in J r that is only 
about half that spanned by the black particles. 

We explain this as a result of the cluster's self-gravity, as fol- 
lows. Particles that are stripped from the cluster mostly escape 
through the Lagrange points L\ and L2. These two points lie near 
the cluster along the line from the barycentre of the host to the 
cluster (Fig. 8.6 of BT08). Thus, it is the particles' radial velocity 
which initially carries them away from the cluster, and it is from 
this velocity component that the particles pay most of the energetic 
penalty for escaping, with the consequence that the radial velocity 
dispersion of the escaping stars is reduced. This reduction in the 
radial velocity dispersion corresponds to a compression in J r of 
the distribution of escaping particles. Once unbound, the particles 
are carried further away from the cluster along a complex trajec- 
tory that ends up with the now-free particle drifting away from the 
cluster according to the mapping in equation (8). Thus, the final 
sum of the energetic penalty is paid from the difference in action 
J Jo between the particle and the cluster, with the net result 
that the unbound distribution uniformly shrinks. The latter effect 
is minor compared to the compression in J r , because much more 
work is done in becoming unbound than in escaping to infinity once 
already unbound. The complete effect is to generate an unbound 
distribution, that looks like the high AJ wings of the pericentre 
distribution, but is compressed in J r and slightly shrunk. 

Looking again at Fig. 10, we note that in the lower-middle 
panel, which corresponds to the 7th apocentre passage, many parti- 
cles have escaped, and the size of the bound distribution has visibly 
shrunk. We understand this as a consequence of the most energetic 
particles having already escaped the cluster, leaving behind a colder 
core. By the time of the 14th apocentre passage, shown in the lower- 
right panel, almost all the particles have escaped. We note that the 
positions of many of the red particles have remained static between 
the lower-middle and lower-right panels. The action-space distri- 
bution of the unbound particles is therefore frozen in place, con- 
firming that self-gravity is unimportant after a stream has initially 
formed. 

In conclusion, we have qualitatively understood the distortion 
of a cluster in action-space as it passes through pericentre and apoc- 
entre along its orbit. We have found that stars stripped at pericen- 
tre form a distribution that is derived from the pericentre distribu- 
tion of bound stars, but is compressed in J r . We have found that 
the pericentre distribution exhibits a high correlation between J r 
and L, and that the gradient of this correlation in (L, J r ) scales 
as ~ y/ J r /L. We typically assume that our cluster orbit will have 
large L and comparatively smaller J r : this would only be untrue 
for extreme plunging orbits, which are not likely to be relevant to 
the problem in hand. Hence the gradient in {L, J r ) will typically 
be less than unity, and the compression will only act to shrink it 
still further: the stripping mechanism produces an action-space dis- 
tribution that is both flattened and very roughly oriented along L. 

6.2. 1 The effect of changing the cluster model or orbit 

We now investigate the qualitative effects on the action-space distri- 
bution of a disrupted cluster of changing the cluster model param- 



eters, or the cluster orbit. The cluster models used in this section 
are CI to C4, detailed in Table 3. These clusters were created ac- 
cording to the schema of Section 6.1, taking the orbit to be 14 in 
the isochrone potential of Section 4.3.1 for the models C1-C3, and 
taking the orbit to be 15 in the same potential for the model C4. 

The cluster model CI, the evolution of which along the orbit 
14 was detailed in the previous section, is used as our baseline to 
which we compare the distributions of the other clusters. The clus- 
ter model C2 has the same profile parameter, W — 2, as does CI, 
but is 10 times more massive. The result is a cluster that is both 
heavier and proportionately larger while being stripped at the same 
galactocentric radius, r s . 

The cluster model C3 has the same mass as does CI, but is 
considerably more concentrated, with a profile parameter W = 6. 
The cluster has an identical truncation radius ft and velocity scale 
a, but has significantly fewer particles near to f t than CI. The par- 
ticles of C3 are generally more tightly bound to the cluster than are 
those of CI. 

The cluster model C4 has the same mass and profile parameter 
as CI, but is specified for the orbit 15, which has lower L than 14, 
and thus a smaller pericentre radius. The resulting cluster is slightly 
more compact, allowing it to survive to a closer galactocentric ra- 
dius than can C 1 . 

A 10 4 particle realization of each of the models C1-C3 was 
placed at a point shortly after apocentre on the orbit 14, and evolved 
forward in time by the FVFPS tree code, using a time step of dt = 
tdyn/100 and a softening length e as specified in Table 3. The total 
period of the simulation was 2.36 Gyr, or almost 7 complete radial 
orbits. Additionally, a 10 4 particle realization of C4 was placed at 
a point shortly after apocentre on the orbit 15, and evolved forward 
in time by the FVFPS tree code, for a total period of 2.21 Gyr, or 
almost 7 complete radial orbits. 

Fig. 12 shows the action-space distribution for the cluster at 
the seventh apocentre passage, for each of these simulations. The 
top-left panel shows C 1 on the orbit 14, and is identical to the lower- 
middle panel in Fig. 10. The top-right panel shows the cluster C2 
on the same orbit; the bottom-left panel shows the cluster C3 on the 
same orbit; and the bottom-right panel shows the cluster CI on the 
orbit 15. 

In the top-right panel, we see that the action-space distribution 
of the more massive cluster is qualitatively identical to that of CI, 
except that the distribution has approximately twice the scale. We 
conclude that cluster mass merely sets the scale of distribution in 
action-space of stripped stars. 

The bottom-left panel shows the distribution from the more 
centrally concentrated cluster C3. In this case, the region occupied 
by bound particles is longer than for CI, reflecting the greater in- 
ternal random velocities of the more concentrated model. However, 
the gross structure of the unbound distribution is approximately the 
same as for CI. This is because the extent of this distribution is 
defined by the least bound stars, stripped from the outer parts of 
each cluster: the behaviour of these stars is determined only by the 
mass of the cluster, and not by its internal configuration. Thus we 
conclude that cluster concentration does not alter the general shape 
of the eventual action-space distribution. 

The bottom-right panel shows the distribution from the cluster 
C4 on the orbit 15, which has the same apocentre radius as 14 but a 
pericentre radius about 33 per cent smaller. Unlike in the other pan- 
els, the cluster has become completely unbound by the 7th apocen- 
tre passage on 15, which is likely to be a result of (r s — r p ) being 
slightly larger for C4 on 15, when compared to the other clusters on 
14, resulting in more efficient stripping at pericentre. 
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Figure 12. The action-space distribution of particles for cluster models on various orbits. Each plot shows the distribution at the 7th apocentre passage. Black 
points are bound to the cluster, red points orbit free in the host potential. From left-to-right and top-to-bottom, the panels show: the CI cluster on the orbit 14; 
the C2 cluster on the orbit 14; the C3 cluster on the orbit 14; and the C4 cluster on the orbit 15. Note the change of scale between successive plots. 



The distribution shown in this plot has approximately the same 
scale in AL as does the distribution from 14, but has approxi- 
mately twice the scale in AJ r . This is a consequence of the gra- 
dient dJ r /dL given by equation (53) being steeper for 15 than 
for 14. Further, since we have already noted that this cluster was 
stripped faster than was CI, the energetic penalty for escaping must 
be lower, and so we expect less compression in J r . The resulting 
distribution is similar to that of the top-left panel, but less com- 
pressed in the J r direction. Thus we conclude that changing the 
cluster orbit can distort the shape of the action-space distribution 
of unbound particles, but does not affect its basic structure. 

6.3 Predicting the stream from the action-space distribution 

We have determined the action-space distribution for several dis- 
rupted cluster models by means of N-body simulation. We now ask 
whether we can accurately predict the real-space path of the stream, 
given the action-space distribution of one of those models. 

Suppose that we know the time since a cluster's first peri- 
centre passage. The angle-space distribution is predicted by equa- 
tion (8), where D is evaluated on the progenitor orbit Jo and 
t — ¥ ta- We will use as an example the action-space distribution 
shown in the bottom-right panel of Fig. 10, which corresponds to 
the 14th pericentre passage of the simulated cluster CI on the orbit 
14. The upper panel of Fig. 13 shows the angle-space configura- 
tion corresponding to this panel: the grey particles are for angles 
directly computed from the results of the N-body simulation, while 
the red particles are for those predicted by equation (8), assuming 
td is known. Also plotted is the frequency vector fio, shown as a 
black arrowed line. The distributions of black and grey particles in 



this panel agree perfectly. Furthermore, both distributions are obvi- 
ously misaligned with the progenitor orbit. 

The lower panel of Fig. 13 shows the real-space configura- 
tion of particles for the same scenario. The grey particles are again 
plotted directly from the results of the N-body simulation, while 
the red particles result from the mapping into real-space of the red 
particles from the upper panel. As in angle-space, the two distri- 
butions agree perfectly. Furthermore, the real-space manifestation 
of the misalignment of the stream with the orbit can be seen: the 
stream delineates a track that has substantially lower curvature than 
the orbit. 

Our attempt to predict the real-space stream configuration 
from the action-space distribution has been completely successful. 
However, any complete model of the bottom-right panel of Fig. 10 
must necessarily be rather complicated. It might not be possible 
to guess the form of this distribution without full N-body mod- 
elling. The dashed line in the bottom-right panel of Fig. 10 is a 
least-squares fit of the action-space distribution to a line. This rep- 
resents a simpler model of the action-space distribution, which it 
might well be possible to guess ab initio for a cluster on a given 
orbit. 

How good a prediction for the stream track can we get from 
this line? The blue lines in Fig. 13 show the results of mapping 
this line into both angle-space and then real-space. It is clearly an 
excellent fit to the stream, in marked contrast to the orbit, which is a 
much poorer model of the stream. Thus, even a very simple model 
of the cluster in action-space — albeit one deduced from an accurate 
knowledge of the distribution — allows us to predict stream tracks 
accurately. 

Finally, we note that during the mapping of this line into real- 
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Figure 13. The distribution of particles for cluster CI, near its 14th apocen- 
tre passage on orbit 14. The upper panel shows the angle-space distribution, 
while the lower panel shows the configuration in real space. The grey par- 
ticles show positions directly computed from the N-body simulation, while 
the red particles show those positions predicted from mapping the action- 
space distribution in the bottom-right panel of Fig. 10. The two distributions 
almost precisely overlap. In both plots, the black arrowed curves show the 
trajectory of the progenitor orbit, while the blue curves show the mapping 
of the dashed line from Fig. 10. The blue curve is everywhere a much closer 
match to the stream particles than is the progenitor orbit. 



space, we need to make a correction to the cluster's action, as de- 
scribed by equation (32), to account for the variation in action down 
the stream. Evaluating equation (46) and equation (49) for the or- 
bit 14, and taking AL ~ 25 kpc kms -1 and dj r ~ 10 kpckms - 
from the bottom-right panel of Fig. 10, we predict errors in the real- 
space of up to ~ 0.2 kpc on account of the finite L distribution, and 
up to ~ 0.15 kpc on account of the finite J r distribution. These er- 
rors would be serious enough to be seen in Fig. 13, and thus the 
correction is required. We note, however, that even these substan- 
tial errors are insignificant compared to the several-kpc discrepancy 
between the stream and the orbit. 



7 THE CONSEQUENCES OF FITTING ORBITS 

We confirmed in Section 4.3 that streams formed in the isochrone 
potential are not well delineated by orbits, and we provided a real- 



istic example where this is the case. Many practical techniques at- 
tempt to use tidal streams to constrain the parameters of the Galac- 
tic potential, by erroneously assuming that such streams do delin- 
eate orbits. In this section, we briefly examine the consequences for 
such techniques of making this faulty assumption. 

Our experiment is based on the simulated tidal stream of 
Fig. 13. This particular example is for a stream at apocentre: this 
is relevant, because many actual observed streams are discovered 
close to apsis, for example the Orphan stream (Belokurov et al. 
2007) and GD-1 (Grillmair & Dionatos 2006). 

We first create two sets of pseudo-data: one corresponding to 
the progenitor orbit, marked as a black line in Fig. 13, and one 
corresponding to the predicted stream track, marked as a blue line. 
Each pseudo-data set contains approximately 40 phase-space coor- 
dinates, sampled evenly from the length of the corresponding track, 
as shown in Fig. 13. 

We now wish to measure the quality with which an orbit for a 
given set of isochrone-potential parameters can be made to fit the 
data. For the purposes of this exercise, we have assumed the func- 
tional form of the potential to be known. Further, we have granted 
ourselves pseudo-data with full and accurate positional and ve- 
locity information. Granting ourselves an unrealistically complete 
pseudo-data set simplifies considerably the matter of finding an or- 
bit that is close to the best fitting one. In practice, one typically 
proceeds with one or more phase-space coordinates unknown, or 
known with substantially less precision than the others. In this case, 
our purpose is solely to demonstrate the errors that can be made by 
naively assuming that an observed stream can be fit with an orbit. 
Granting ourselves an unrealistically complete set of pseudo-data 
does not diminish our conclusions in this regard. 

For each set of potential parameters, we choose an orbit as 
follows. We first select a datum near the centroid of the stream, 
and declare that our chosen orbit must pass directly through this 
datum. Although it may be that some nearby orbit, one that does 
not pass directly through this point, would make a better-fitting or- 
bit, any such orbit must pass very close to the selected datum, be- 
cause it is close to the centroid. Thus, such an orbit cannot be much 
better-fitting than one that passes directly through the datum. Hav- 
ing chosen a datum, for a given set of potential parameters, an orbit 
is defined. 

Having chosen our orbit, a goodness-of-fit statistic x 2 is cal- 
culated as follows. For each datum in the stream, with phase-space 
coordinate w„ a location along the orbit w< is chosen that mini- 
mizes the square difference 



(W. - Wj) . 

Having obtained the the goodness-of-fit \ 2 is defined by 



(54) 



(55) 



where j are the phase-space coordinates, and Oj is the rms of 
(wij + w'ij)/2 over i. This \ 2 statistic provides a dimensionless 
measure of the phase-space distance between the best-fitting orbit 
in a given potential, and the pseudo-data. 

If the pseudo-data set were a sample of a perfect orbit in some 
potential, we expect the value of x 2 to be exactly zero, when the 
correct potential parameters are used. As the potential parameters 
are varied away from their true values, we expect the value of x 2 
to rise, as the best-fitting orbit becomes a steadily worse represen- 
tation of the data. Hence, we expect minima in x 2 to be associated 
with the potential parameters that are optimum, from the perspec- 
tive of fitting an orbit to the data. We seek such minima by plotting 



16 A. Eyre and J. Binney 




10 20 30 40 50 

GM/(207kmj-') 2 kpc 

260 F H 




5 10 15 20 

r / kpc 



Figure 14. Top panel: goodness-of-fit \ 2 f° r eacn pseudo-data set, to the 
best-fitting orbit in the potential. The colours of the curves for each pseudo- 
data set identify it with the corresponding curve from Fig. 13 from which it 
was derived. Bottom panel: rotation curves for the optimum potential found 
for each pseudo-data set. In all cases, we have required that the circular 
velocity v c = 240 kms -1 at i?o = 8kpc. This reduces the potential 
to one of a single parameter; in this case, mass. The quality of the fit to 
the blue data set is significantly degraded, compared to that of the black. 
Utilising the high-curvature, blue data set as a proxy for an orbit causes us 
to overestimate the host mass by approximately 21 per cent. 

the value of x' 2 over a range of likely values for the potential pa- 
rameters. 

7.1 Results 

The upper panel of Fig. 14 shows the goodness-of-fit \ 2 > f° r the 
best-fitting orbit, in an isochrone potential with mass parameter 
GM as shown. For each value for GM, the isochrone parameter 
b is chosen such that the potential reproduces the fiducial rotational 
speed of v c — 240 kms -1 at Rq = 8 kpc, since in practical usage 
one would generally require any acceptable potential to reproduce 
other observed features of the Milky Way galaxy, such as the circu- 
lar velocity at the Solar radius. This has the effect of reducing the 
dimensionality of the search for the best-fitting solution to a single 
parameter. 

Along each of the black and blue curves in Fig. 14, x 2 has a 
minimum at the potential parameters for which the best-fitting orbit 
has been found. In the case of the black curve, which attempts to fit 
to pseudo-data derived directly from the progenitor orbit, the value 
of GM is correctly identified with high precision, thus validating 
the technique. 

For the blue curve, which attempts to fit pseudo-data derived 
from the stream track, the quality of the fit at optimum GM is sig- 
nificantly degraded when compared to the black curve. This indi- 
cates that, although an optimal value for the mass parameter is be- 



ing found, the fit to the orbit is still not perfect there. Furthermore, 
the value of GM for which the optimal fit is found is 21 per cent 
larger than the true value. 

The attempt to constrain the potential parameters by using 
misaligned streams as proxies for orbits has led us to err. The high- 
curvature of the stream, relative to the orbits of the stars that com- 
prise it, has incorrectly led the fitting procedure to a significantly 
larger enclosed mass, in order to provide the larger gravitational 
force necessary to produce a more highly-curved orbit. This can 
clearly be seen from the lower panel of Fig. 14, which shows the ro- 
tation curves for the optimal potentials for both sets of pseudo-data. 
In the case of the stream-derived pseudo-data, the circular speed of 
the optimum potential is significantly higher than the truth at those 
distances where our stream segment lies. 

In conclusion, we find that there is a risk of substantial system- 
atic errors in parameter estimation being made, if one attempts to 
constrain the potential using streams, and one assumes that streams 
perfectly delineate orbits. In practice, one would hope that the sig- 
nificant degradation of the minimum \ 2 °f the fit to the stream- 
derived data, when compared to orbit-derived pseudo-data, would 
be noticeable. However, it is often the case that random noise in 
the measurements are sufficient to reduce the quality of fits such 
that such discrimination is not possible (Binney 2008; Eyre & Bin- 
ney 2009b). If this is so, then even with the perfect error-free data 
used in this section, the uncertainty on the deduced mass parameter 
could not be better than ±20 per cent. In practice, the uncertainty 
is likely to be much, much worse. Only by fitting stream data to 
physically-motivated stream tracks instead of orbits can we hope 
to overcome this limitation, and obtain the tantalising diagnostic 
precision proffered by the analysis of streams in such as Binney 
(2008). 



8 NON-SPHERICAL SYSTEMS 

We have investigated the formation of streams in spherical poten- 
tials, and demonstrated that they are not necessarily well modelled 
by orbits, but that with some prior knowledge of the system, we can 
accurately predict their tracks. However, many real stellar systems 
in the Universe are not spherical. In particular, our own Galaxy, 
whose potential we are interested in probing with streams, is sig- 
nificantly flattened. 

In the case of spherical systems, the stream can only be mis- 
aligned with the progenitor orbit within the plane to which both 
stream and orbit are confined. Orbits in flattened potentials are 
generically three-dimensional as a result of the instantaneous or- 
bital plane precessing around the potential's symmetry axis. Differ- 
ent parts of the stream have orbit planes which precess at different 
rates, so the stream is not strictly confined to a plane. Moreover, the 
plane that most nearly contains the stream generally differs from 
the instantaneous orbit plane of the progenitor. These complexi- 
ties come on top of the misalignments within the best-fitting plane 
that we have studied in spherical potentials. They could well be the 
reason that all attempts to find a single orbit that fits the Sagittar- 
ius dwarf stream have had limited success (Fellhauer et al. 2006; 
Niederste-Ostholt et al. 2010). In this section, we investigate the 
formation of streams in a flattened Stackel potential. 

8.1 Stackel potentials 

Regular orbits in non-spherical potentials can be described by a sys- 
tem of three actions, but the actions and angles can be computed by 
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standard means only if the potential is a Stackel potential (§3.5.3 
of BT08). An exhaustive treatment of these potentials is given in 
de Zeeuw (1985) and de Zeeuw et al. (1986), from which we adopt 
our notation. The example we consider is an oblate axisymmetric 
potential, for which the appropriate coordinate system is that of 
prolate spheroidal coordinates, (A, <f>, u); A is constant on confocal 
ellipsoidal surfaces of revolution, while v is constant on the hy- 
perbolic surface of revolution that cut the ellipsoids orthogonally. 
Hence A is a radial coordinate and v is analogous to latitude. The 
coordinate system is specified by two negative quantities (0,7). 
The potential has the form 

(A + 7 )G(A)-(, + 7 )GM 



\ — V 



(56) 



where de Zeeuw's function G(r) is determined once the density 
profile p(z) along the z-axis has been chosen (see e.g. equation 
(23), de Zeeuw et al. 1986). Thus, the model is completely specified 
by p{z) and the scaling parameters (a, 7). The actions for Stackel 
potentials of this type are ( J\, L z , J„), corresponding to motion in 
(\,<t>,v) respectively. 

8.2 Galaxy models with Stackel potentials 

de Zeeuw et al. (1986) shows that if one requires the density every- 
where to be non-negative, it is not possible to formulate a Stackel 
model in which the density p(r) falls off with distance from the 
z-axis more rapidly than r~ 4 as r — > 00. This is because an ele- 
mentary density on the z-axis p(z) = S(z — zo) requires an off- 
axis density term that falls as r~ 4 . This behaviour rules out many 
classes of galaxy models, including discs with exponential density 
profiles. However, we can construct models in which the density 
falls more slowly than r -4 as r — > 00. In particular, models with 
asymptotically flat rotation curves, i.e. those in which p(r) ~ r~ 2 , 
are allowed. 

In the models used in this section, we specify the z-axis den- 
sity profile 



Pz(z) 



"7Po 



-7Po 



(z2- 7 ) 



(57) 



where we have made use of z 2 = r + 7. In this case, de Zeeuw's 
function G(r) can be written in closed form (equation (49) with 
s = 2, de Zeeuw et al. 1986). Models specified by equation (57) 
become spherical at large radii and have a rotation curve that is 
asymptotically flat, with 



lim v c = — 4/7rGpo7- 



(58) 



In the core of these models, the surfaces of constant density are 
approximately ellipsoidal, with axis ratio 1 



2q 2 



(i- q *y 



(1 



2 , 

q + 



2 , 2\ 

2 log q J 



(59) 



where the central potential axis ratio q — y/j/a. The models are 
completely specified by choosing a shape via q, a mass scale via po, 
and a distance scale via 7. The combination of an asymptotic loga- 
rithmic 'halo' and a flattened 'disc' in these models allows them to 
make a fair representation of the observed properties of disc galax- 
ies, although the lack of freedom in the models severely restricts 
the shape of the density distribution that can be achieved. 

1 Equation (46) in de Zeeuw et al. (1986) presents a formula for this quan- 
tity, but it is incorrect. 
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Figure 15. The upper panel shows contours of log p/po for the Stackel 
model SP1, utilized in this section. The lower panel shows contours of the 
potential <E> for the same. 



The potential used in this section, SP1, was chosen to simulate 
a galaxy with a massive disc. The axis ratio of the density distribu- 
tion is fixed to be 10 near the solar radius, 7?o = 8kpc, which is 
approximately the same ratio as for the (exponential) thin disc pro- 
file of the Milky Way (see Table 2.3 of BT08). The specification is 
completed by requiring the rotation curve to peak at _Ro = 8 kpc 
with a circular speed v c = 240 km s -1 . The resulting parame- 
ters of the SP1 potential are: p = 3.61 x lO 9 M kpc _3 ,a = 
-29.64kpc 2 , 7 = -8.89 x lO^kpc 2 . 

We note that the asymptotic circular velocity in this model is 
v c — 42kms _1 , which can be regarded as the halo contribution 
to the circular speed. The model is too centrally concentrated, and 
the halo contribution is too weak, to realistically model the Milky 
Way. However, it is highly flattened, and so makes an interesting 
example in which to study stream geometry. Contours of constant 
density and potential for this model are shown in Fig. 15. 



8.3 Stream misalignment in the Stackel potential SP1 

We now consider the geometry of streams formed in the Stackel 
potential SP1. Although the actions J can be defined in terms of an 
integral over a single coordinate, closed-form expressions for J do 
not exist. Instead, the integrals in the expressions for J have to be 
evaluated numerically. Similarly, expressions for both the frequen- 
cies ft and their derivatives Vjfi can be derived, but not in closed 
form, and the integrals that they contain must also be evaluated nu- 
merically. 

The expressions for J and CI appear in de Zeeuw (1985). The 
detailed expressions for Vjf2 are algebraically rather involved, and 
are presented in the appendix of Eyre (2010b). Here, we simply 
note that having evaluated these quantities for a particular orbit Jo, 
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Figure 16. Details of the stream geometry for the Stackel potential SP1. Left panels: contours for the misalignment angle i9, in degrees, between the principal 
eigenvector of D and fin, shown as a function of J. Right panels: contours for the eigenvalue ratio A1/A2. The top panels show the plane in action-space 
with J v = 20.7 kpc km s -1 , while the bottom panels show the plane in action-space with L z = 414 kpckms -1 . The range of actions covers a variety of 
interesting orbits: details of orbits at the extremes of the range are given in Table 4 



Table 4. The coordinate extrema of selected orbits from Fig. 16, illustrating 
the variety of orbits covered by that figure. The actions are expressed in 
kpc km s~ 1 , while the apses are in kpc. 
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the eigenvectors e„ and the eigenvalues A n are computed directly 
from the matrix D(Jrj) = Vjf2| Jo by standard methods. 

The left panels of Fig. 16 show contour plots of the misalign- 
ment t? in angle-space between the principal direction of D and the 
frequency vector firj, where # is calculated from equation (27), as 
was the case for systems of two actions. The right panels of the 
same figure show contours of the ratio A1/A2. The range of actions 
shown in these plots covers a variety of interesting orbits; the apses 
of the orbits at the extremes of the range are described in Table 4. 

As with the equivalent plots for the isochrone potential 
(Fig. 3), we see that the principal direction of D is never perfectly 
aligned with firj. In this very flattened potential, streams with low 
J„ have the greatest degree of misalignment, at about ~ 15 deg. 
These orbits never go far from the disc. The misalignment dimin- 



ishes with increasing J„, falling to ~ 4 deg for orbits with apses 
in z of some tens of kpc. Hence, in this very flattened poten- 
tial, there is much more prospect for dramatic angle-space mis- 
alignment than the small misalignments we encountered with the 
isochrone potential (Fig. 3). Fig. 16 also shows that the eigenvalue 
ratio A1/A2 > 10 everywhere for the SP1 potential; thus, we con- 
clude that highly elongated streams will form on all orbits which 
permit a cluster to be disrupted. 

As with the spherical case, the precise behaviour of any given 
stream depends on both the potential and the action-space distri- 
bution of its stars. To proceed further we must again consider a 
specific example, by means of N-body simulation. 

8.4 A stream in the Stackel potential SP1 

Fig. 17 shows the real-space trajectory of the orbit SOI, evaluated 
in the Stackel potential SP1. This orbit has actions (J\, L Z ,J V ) = 
(252.3, 2618., 20.4) kpckms" 1 and apses of approximately R = 
(8, 18) kpc in the galactic plane, and z — (—2,2) above and below 
the plane. It is thus fairly representative of an eccentric orbit that 
might be occupied by a globular cluster orbiting close to a galactic 
disc. 

The cluster model C5 (Table 3) is a King model specified for 
the orbit SOI according to the schema of Section 6.1. The model 
has the same mass and profile parameter as CI, and is very similar 
in all other attributes, because the orbit SOI is similar to the orbit 
14 for which CI was specified. A 10 4 particle realization of the 
C5 was made by random sampling of the King model distribution 
function. This cluster was placed close to apocentre on the orbit 
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Figure 17. The real-space trajectory of the SOI orbit in the Stackel potential 
SP1. The upper panel shows a plan view of the galactic plane, while the 
lower panel plots height above the plane, z, against azimuthal coordinate, 



SOI and evolved forward in time by the FVFPS tree code, with 
time step At — tdyn/100 and softening parameter e as specified 
in Table 3. The total period of the simulation was 2.15 Gyr, or 7 
complete radial oscillations. 



8.4.1 Action-space distribution 

Fig. 18 shows the evolution in time of the action-space distribution 
of this simulated cluster. Each column of panels shows the distribu- 
tion at a different point in time. The upper panel in each row shows 
the orthographic projection of the actions onto the (L z , J\) plane, 
while the lower panel of each row shows a similar projection onto 
the (L z , J„) plane. In all panels, the appropriate projection of the 
mapped frequency vector, D _1 fio, is shown as an arrowed black 
line. 

The left hand column shows the actions when the cluster is 
near to its first pericentre passage. In the upper panel, the distribu- 
tion is somewhat flattened, and oriented with positive gradient in 
AJ\/ AL Z . This behaviour is analogous with that seen in the top- 
right panel of Fig. 10: the motion of the cluster is predominantly in 
(A, (f>), thus ( J\, L z ) are good proxies for the radial action J r and 
angular momentum L respectively. J\ and L z are therefore highly 
correlated for a cluster near apsis on this orbit, in analogy with the 
mechanism described by equation (53) in Section 6.2. 

Meanwhile, the distribution in the lower panel is both narrow 
and oriented almost exactly along L z . We can understand the shape 
as follows. For this orbit, which is confined to be close to the plane, 
Jv ~ Jz/2, where the factor of 2 appears because J„ is defined on 
a path restricted to only one side of the plane. J z can be estimated 
by close analogy with equation (29). Hence, the spread in J„ for a 
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where R p is the galactocentric pericentre radius in cylindrical 
coordinates. Evaluating this expression for the orbit SOI gives 
AJ„/AL Z ~ 0.04, which we see from the lower-left panel of 
Fig. 1 8 is close to exact. 

The flat orientation of the lower-left panel we explain by 
pointing out that, as Fig. 17 shows, the motion in v in this example 
is almost decoupled from the radial motion. This means that the v 
coordinate need not be at apsis when the cluster is at pericentre, 
and thus the arguments of Section 6.2, which force a correlation 
between J\ ~ J r and L z ~ L near pericentre, do not apply. For 
an orbit in which J„ is more strongly coupled to the radial motion, 
we would expect to see the characteristic tilting of the (L z , J„) dis- 
tribution near pericentre and apocentre, as a correlation between J„ 
and L is forced. 

The middle column of panels in Fig. 18 shows the action- 
space distribution at the subsequent apocentre passage. Bound and 
unbound particles are shown in black and red, respectively. The 
action-space structure of the unbound particles in the upper panel 
bears striking similarity to that shown in Fig. 10, as might be ex- 
pected since (J\,L Z ) make good proxies for (J r ,L). The same 
physical principle for the disruption of the cluster applies here as 
in the spherical case; that is, the particles will escape the cluster 
through the Lagrange points L\ and L2 by first travelling radially, 
so the action-space distribution will be compressed in this direc- 
tion. Thus, the range of AL Z and AJ V for the unbound stars is 
about the same as in the left hand panels, but the range of A J\ is 
markedly less. 

The right hand panels of Fig. 18 show the action-space distri- 
bution at the 7th apocentre passage. The structure is essentially the 
same as that of the middle panels, except that all particles are now 
unbound. Also plotted in each of the right hand panels is a dashed 
line, which has been least-squares fitted to the unbound particles. 
We note that the image of the frequency vector and the dashed line 
are highly misaligned in the upper-right plot; because of this, we 
expect the stream to be significantly misaligned with flo in angle- 
space. 

In conclusion, we find that in very flattened potentials, dis- 
rupted clusters form an action-space distribution that is wholly 
analogous with, although necessarily more complicated than, that 
found for disrupted clusters in spherical potentials. 



8.4.2 The effects of disc shocks 

Unlike with a spherical potential, an axisymmetric potential allows 
for tidal forces other than those felt during pericentre passage to 
act upon an orbiting cluster. In particular, the passage of a cluster 
through a massive galactic disc will subject a cluster to a tidal force 
that is of comparable magnitude to that felt when close to pericen- 
tre. 

The tidal stress imposed on a cluster at pericentre has a tensile 
component, which acts to strip stars from the cluster. Conversely, 
the tidal stress imposed by a disc passage is entirely compressive in 
nature. Hence, stars are not actively stripped from a cluster during a 
disc passage. Instead, the action of such 'disc shocks' is to heat the 
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Figure 18. Action-space distribution for the simulated cluster C5, on the orbit SOI in the Stackel potential SP1, at various points in time. The upper panels show 
the orthographic projection of distribution onto the (L z , JpJ plane, while the lower panels show a similar projection onto the (L z , J„) plane. Distributions 
are shown at the following times: left panels, the first pericentre passage; middle panels, the subsequent apocentre passage; right panels, the 7th apocentre 
passage. Black particles are bound to the cluster, while red particles orbit free in the host potential. The dashed lines in the bottom panels are lines that have 
been least-squares fitted to the particles. 



cluster, perhaps repopulating the outer edges of the cluster, which 
were depopulated during a previous pericentre encounter (Spitzer 
1987, §5.2a). 

The net effect of disc shocks on the stripping process is a 
faster and more complete disruption of the cluster than would take 
place for an unshocked cluster exposed to equivalent pericentric 
tidal stress. Since the vast majority of stars continue to be stripped 
at pericentre even when the effect of disc shocks is significant, 
the gross action-space distribution resulting from the stripping of 
a shocked cluster will remain as previously described. However, 
since the disc shocks act to increase the velocity dispersion of the 
cluster between stripping events, it is likely that the wings of the re- 
sulting action-space distribution will be populated with more stars 
than would otherwise have been the case. 

8.4.3 Predicting the stream track 

Fig. 19 shows the angle-space configuration for the simulated clus- 
ter near its 7th apocentre passage. The grey particles are for angles 
that have been computed directly from the output of the N-body 
simulation. The arrowed black line is fin, while the blue line is the 
map of the dashed-line from Fig. 18. In both projections, the blue 
line is a much superior match to the data than is the orbit. 

Fig. 20 shows equivalent plots to those of Fig. 19, but in real 
space. The misalignment between the stream and the progenitor 
orbit in angle-space is seen to map into a large misalignment in real 
space. Attempting to constrain halo parameters by fitting orbits to 
the stream shown in Fig. 20 would not produce sensible results. 



Conversely, the map of the dashed-line model for the action-space 
distribution, shown in Fig. 18, clearly provides an excellent proxy 
for the track of the stream. We note that, as per our expectation, 
there is significant distortion both in and out of the effective orbital 
'plane' . 



9 CONCLUSIONS 

We have studied the mechanics of the disruption of low-mass clus- 
ters, and the subsequent formation of tidal streams. In particular, 
we have examined the conditions required for the formation of tidal 
streams from such clusters, and to what extent such streams delin- 
eate the orbits of their stars. With regard to the latter, we are moti- 
vated by techniques that promise to constrain the form and param- 
eters of the Galactic potential when observations of thin streams 
are analyzed under the assumption that such streams precisely de- 
lineate orbits (Binney 2008; Eyre & Binney 2009a; Willett et al. 
2009; Newberg et al. 2010; Koposov et al. 2010). 

We utilize action-angle variables extensively in our approach 
to the problem. These coordinates allow a convenient and natural 
description of the physical processes that occur in cluster disruption 
and stream formation. The stream has a structure in action-space 
that is formed at the point of disruption, and subsequently frozen 
for all time. The angle-space structure of the stream is the image 
of that action-space distribution under a linear map, which is itself 
a function of the host potential and the progenitor orbit only. This 
angle-space image elongates linearly with time, corresponding to 
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Figure 19. Angles for the N-body cluster shown in the right hand panels 
of Fig. 18. The blue line shows the predicted stream, resulting from the 
mapping of the dashed line from the same plots. The blue line is clearly 
a much better representation of the stream than is flo , represented by an 
arrowed black line. 
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Figure 20. Real-space representation of an N-body simulation of cluster 
model C5 evolved on the orbit SOI, in the SP1 potential. The points are 
shown near the 7th apocentre passage following release. The black line 
shows the trajectory of an individual orbit; the blue line shows the predicted 
stream path. 



the extension of the stream, but its shape remains otherwise un- 
changed. Disrupted clusters elongate to form streams because this 
linear map preferentially stretches the angle-space structure along 
a particular direction. We find this preferential stretch to occur for 
likely orbits in all the realistic potentials that we consider: the only 
exception is the spherical harmonic oscillator potential, in which 
clusters cannot be stripped at all. 

The real-space structure of a stream is a function of the cor- 
responding action-angle structure and the host potential. We have 
investigated the conditions required for a stream to be perfectly de- 
lineated in real space by its progenitor's orbit. One necessary con- 
dition is that the angle-space stream be perfectly delineated by the 
frequency vector. 

We have shown that this will always be the case for a stream 
formed in a Kepler potential. Conversely, we have shown that this 
is generally not the case for a stream formed in an isochrone po- 
tential: angle-space streams will generally be misaligned with the 
frequency vector by a few degrees. Since the isochrone potential 
is a more realistic representation of a galactic potential than is the 
Kepler potential, we infer that streams formed in galactic potentials 
will not generally be perfectly delineated by orbits. Indeed, it would 
appear that the alignment of streams with orbits is a property that 
is peculiar to the Kepler potential. Since streams exist in the outer 
parts of galaxies where the potential has a substantial monopole 
component, we propose this as an explanation for the observation 
that streams appear to be approximately delineated by the orbits of 
their stars. 

We have examined the real-space manifestation of angle-space 
streams that are misaligned with their frequency vector. We find 
that the real-space effect depends on the orbital phase at which the 



stream is observed: at apsis, the stream track will have a curvature 
that differs from that of the progenitor orbit, while away from apsis 
the stream and its orbits will be spatially misaligned. 

We have also examined the real-space effects of a finite action- 
space distribution for a stream. We conclude that, with a cluster de- 
scribed by a large, but realistic, velocity dispersion, knowledge of 
the angle-space structure and the progenitor orbit is no longer suf- 
ficient to accurately predict the real-space track of a stream. We de- 
scribe a first-order correction to the actions of the stream stars that 
is deduced from observables, plus an estimate of the time elapsed 
since the initial cluster disruption event. This correction suffices 
to allow accurate predictions of real-space tracks to be made for 
streams formed from larger clusters. 

We have confirmed the validity of our analyses with the use 
of N-body simulations of the disruption of clusters in the isochrone 
potential. These show that disrupted clusters indeed form streams 
that are poorly represented by the trajectory of the progenitor orbit. 

We have used these simulations to inform an explanation of 
the action-space structure of disrupting clusters. We find that the 
action-space distribution of a disrupted cluster takes a character- 
istic shape, the details of which depend on the cluster model, the 
cluster orbit and the host potential. We are able to explain all fea- 
tures of this distribution in terms of the basic physical processes 
that apply to clusters undergoing tidal disruption. We further show 
that by utilising a simple, straight-line model of the action-space 
distribution, we are able to predict the real-space stream track of a 
stream with very high accuracy. 

We have briefly examined the consequences of our findings for 
techniques that attempt to constrain the Galactic potential under the 
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assumption that streams can be fit with orbits. We utilize such an 
algorithm to attempt to constrain the mass of an isochrone poten- 
tial, subject to the rotation speed at the Solar circle being equal to 
a fiducial value, and subject to input pseudo-data being fit with an 
orbit. We validate our optimisation procedure using input pseudo- 
data for a track that perfectly delineates an orbit. When the input 
pseudo-data are for a simulated stream in which the track shows 
an increased curvature with respect to orbits in the correct poten- 
tial, the algorithm fails to identify the correct parameters: in our 
example, the response of the algorithm was to report an optimum 
mass parameter 21 per cent greater than the true value, albeit with 
a substantially decreased quality-of-fit compared to attempts to fit 
to pseudo-data representative of a true orbit. We conclude that large 
systematic errors can be made when attempting to optimize poten- 
tial parameters under the assumption that streams act as proxies for 
orbits. 

We have examined stream formation in a flattened, axisym- 
metric Stackel model, which provides a fair approximation to the 
potential near a massive galactic disc. We find that in this poten- 
tial the principal direction of the linear map is misaligned with the 
frequency vector by ~ 10 dcg. Hence in highly flattened potentials 
like this, streams will not be well represented by orbits. 

We performed an N-body simulation of a disrupting cluster on 
an approximately planar orbit in the flattened potential. We found 
that the angle-space misalignment between the stream and the fre- 
quency vector is indeed large, and that the resulting real-space track 
is very poorly represented by the progenitor orbit. However, we 
again find that a simple straight-line model of the action-space dis- 
tribution predicts the corresponding real-space stream track with 
superb accuracy. We also find that in our Stackel example, the 
action-space distribution resulting from the disruption of clusters 
is directly comparable to that found in the isochrone potential, con- 
firming the generality of our earlier conclusions. 

The consequences of our findings for potential-optimization 
techniques are significant: large systematic errors are possible if 
orbits are fitted to streams. Before attempting such activity, it is 
now necessary to check whether the segment of the stream in ques- 
tion could be legitimately represented by an orbit. Similarly, is also 
important for this check to be performed retrospectively for those 
streams — e.g. Orphan stream analysis of Newberg et al. (2010) and 
the GD-1 analysis of Koposov et al. (2010) — that have already been 
used to constrain the parameters of the Galactic potential, since the 
validity of those analyses is now in question. 

Such a check could always be performed using N-body simu- 
lation, where it is believed that the progenitor model, the orbit and 
the applicable potential are well known. However, the application 
of the techniques in this paper would constitute a more general ap- 
proach, in which the implications of an incorrect model or orbit can 
be readily addressed. 

The techniques presented in this paper depend on the ability 
to compute action-angle variables from conventional phase-space 
coordinates, and vice versa. In this paper, we have restricted our 
consideration to models for which global transformation relations 
can be written down. The only such models are those with poten- 
tials of the Stackel type, which includes all spherical potentials as a 
limiting case (de Zeeuw 1985). However, it is difficult to construct 
a realistic flattened galaxy model that is satisfactory in all respects 
whilst restricted to using only Stackel potentials. In general, one 
would like to work with more sophisticated models for which no 
general transformation between action-angle variables and phase- 
space coordinates can be made. 

Fortunately, the "torus machine" (McMillan & Binney 2008) 



enables the actions, the frequencies and their derivatives to be accu- 
rately and quickly computed for regular orbits in realistic Galactic 
potentials. By utilising the torus machine, the techniques explored 
in this paper could be extended to work with such models. This ap- 
proach could be used to investigate the impact on streams of orbital 
resonances and chaos: with the torus machine one can construct 
approximate angle-action variables even for potentials that support 
chaos (Kaasalainen & Binney 1994). Resonances and chaos play 
rather a small role in the dynamics of axisymmetric galactic poten- 
tials, but one might expect the sensitivity of chaotic orbits to their 
initial conditions to lead to a measurably more rapid spread of tidal 
streams in chaotic rather than integrable regions of phase space. 

If a check — however performed — confirms that the stream is 
well modelled by an orbit, then diagnosing the potential by the fit- 
ting of orbits to streams is appropriate. However, if the stream is 
not well modelled by an orbit, as will generally be the case, the 
technique of fitting orbits should not be used. 

In such circumstances, one might resort to N-body shooting 
methods to compute stream tracks to feed to an optimization algo- 
rithm (Johnston et al. 2005). However, the computational expense 
of such a technique would be a major burden, and would severely 
limit the quality of the optimization. 

The results of this paper present a possible alternative. We 
have found that with simple models of the action-space distribution 
of a disrupted cluster, such as can be readily obtained from a single 
N-body simulation, we can reliably and accurately predict the track 
of a stream, even when it diverges significantly from the trajectory 
of the progenitor orbit. In principle, we can compute these tracks 
with no more computational effort than it takes to integrate an orbit. 
Existing techniques for constraining potential parameters by fitting 
orbits could be readily adapted to fit stream tracks instead. 

To achieve this goal, the key hurdle is the development of a 
predictive theory for the action-space structure of a disrupted clus- 
ter, valid for any problem parameters of our choice. A possible ap- 
proach would be to use N-body simulations to obtain the action- 
space distribution for a small number of cluster models on a set 
of possible orbits, for a range of potentials. The resulting distribu- 
tions could be used as a basis set, which would be interpolated and 
distorted to provide an estimate for the action-space structure for 
any given cluster on any chosen orbit. The required distortions for 
changes to cluster model and orbit parameters have already been 
touched upon in this paper, although in order to be of practical use, 
a complete quantitative theory of these distortions will be required. 

In the future, the techniques presented here may well be ap- 
plicable to the Sagittarius dwarf stream. In this paper we have as- 
sumed that our low-mass clusters do not affect the host potential. 
This may not be true in the case of a heavy Sagittarius progenitor. 
Further study of the effect of a live host potential on the mechanics 
of stream formation will be required for the techniques to be reli- 
ably applicable to the Sagittarius dwarf stream. Also, the action- 
space structures discussed will be, in part, composed of dark mat- 
ter in the case of streams formed from dwarf spheroidal galaxies. 
Further investigation of the stripping of satellites using models in- 
cluding both dark particles and stars will be required to determine 
which parts of the stream will remain visible, and which will be 
composed solely of dark matter. 
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